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ABSTRACT 

A development of perturbations of number densities of ions and electrons during 
recombination epoch is analysed. The equations for relative perturbations of ioniza- 
tion fractions were derived from the system of equati ons for accurate computat ion of 
the ionization history of the early Universe given by ISeager et all lIlflflflL 1200(1) . It is 
shown that strong dependence of ionization and recombination rates on the density 
and temperature of plasma provides the significant deviations of amplitudes of ion- 
ization fractions relative perturbations from ones of baryon matter density adiabatic 
perturbations. Such deviations are most prominent for cosmological adiabatic pertur- 
bations of scales larger than sound horizon at recombination epoch. The amplitudes of 
relative perturbations of number densities of electrons and protons at last scattering 
surface exceed by factor of ~5 the amplitude of relative perturbation of baryons total 
number density, for helium ions this ratio reaches the value of ~18. For subhorizon 
cosmological perturbations these ratios appear to be essentially lesser and depend on 
oscillation phase at the moment of decoupling. These perturbations of number den- 
sities of ions and electrons at recombination epoch do not contribute to the intrinsic 
plasma temperature fluctuations but cause the "corrugation" of last scattering surface 
in optical depth, 5zdec/ {zdec + 1) ~ — ^fc/S, at scales larger than sound horizon. It may 
result into noticeable changes of precalculated values of CMB polarization pattern at 
several degrees angular scales. 
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background 



INTRODUCTION 



Cosmic microwave background (CMB) radiation coming from recombination epoch has become one of the most powerful 
observational probes for cosmological models of our Universe and formation of its large-scale structure. Indeed, the full-sky 
maps of cosmic microwave temperature fluctuations obtained by Wilkinson Microwave Anisotropy Probe (WMA P) during flrst 
year of observations have given a possibility to determine the cosmological parameters with high accuracy, ~ 2% tBenn et et alJ 



el 



Mm 



erde et alllinnl : ISnergel et 3,11120031 ) . The current data from three year of WMAP observations 
and, especially, expected from future mission Planck, improve precision of the CMB power spectrum determination to the 
level of accuracy of numer ical prec alculations done by the most advanced codes in the framework of specific model. For 
example, CMBfast code bv iSeliak fc Za ldarriaga (1996); Zaldairriaga_&_Seljak ( 1 99^ ) has intrinsic accuracy ~ 1%. An adequate 
calculation of the recombination process is crucial for modelling the power spectrum of CMB temperature fiuctuations and 
polarization. 
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199Q; Rubicki & Dell'Antonio 1993 and citine therein) the main processes have been studied 



using the 3-level approximation of hydrogen and helium atoms. An accuracy of few percents has been achieved. The most 
complete analysis of cosmological recombination processes with taking into account the multi-level structure of hydrogen and 
helium atoms (~ 300 levels) and non-equilibrium ionization-recombination kinetics has been performed by Seager, Sasselov 
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& Scott (2000). Also all known plasma thermal p rocesses were take n into account therein. These authors have provided 
cosmological community with software RECFAST iSeager et alJllflflflll which ensures an accuracy of calculation of number 
density of electrons ~ 1%. This code was used by number of authors to calculate the transfer function of density perturba- 
tions and power spec t rum of CMB temperature fluctuations and polarization, in particular, by SMiak_& Zaldarriaga (199^; 
IZaldarriaga fc SeliakI fippfll for publicly available software CMBFAST. However, the researches aimed on improving the cal- 
culation of recombination and decoupling of the thermal radiation from baryon pla sma are still going on (see recent pa pers by 
Ipubrovich & Grachev (2005); Chluba & Sunvaev (2006); Kholupe nko. Ivanchik fc Varshalovich (2005); Wong. Scager fc Scotd 
l|2fl05f) and citing therein). 

In this paper the evolution of number density perturbations of hydrogen and helium ions and electrons in the field of 
cosmological matter density perturbations is studied. The cosmological adiabatic perturbations closely associate the density 
and temperature variations in baryon-radiation component, leading to the corresponding variations of photorecombination 
and photoionization rates, which in turn can cause appreciable deviations of relative perturbation amplitudes of number 
density of ions and electrons from corresp onding ampli t ude fo r total number density of baryon nuclei. The effective 3-level 
models of hydrogen and helium atoms by ISeager et al.l l|l999f) and their software RECFAST were used as a basis. Here all 
calculations were carried out for z > 100 when recombination and dissociation processes of hydrogen negative ions H~ and 
molecules H2 and Hj can be neglected due to their insignificance. 

In the first section the basic equations for hydrogen and helium recombination in homogeneous expanding Universe 
are presented along with their numerical solutions for ACDM model. The definitions of relative perturbations of number 
densities of ions and electrons, their properties, equations for their evolution and results of integration for the stationary mass 
density and temperature initial perturbations are given in the second section. The third section is devoted to the analysis of 
perturbations of electron number density within the adiabatic mass density perturbations of different scales in ACDM model. 
An estimation of contribution of this effect into the "corrugation" of last scattering surface and CMB temperature fiuctuations 
is carried out in the 4th section. 



1 COSMOLOGICAL RECOMBINATION OF HYDROGEN AND HELIUM ATOMS: DEFINITIONS, 
EQUATIONS, RESULTS 

Let us introduce the following definitions: nni and nnii denote the number densities of neutral and ionized hydrogen atoms 
respectively, riHoi, ?ihgII and nHciii - number densities of neutral, singly and double ionized helium atoms, no — nHii + nHcii + 
2nHciii - number density of free electrons, jih = WHi+nnii - total number density of hydrogen nuclei, riHo ~ n-Hei+n-Heii+n-Heiii 

- total number density of helium nuclei. It is conveniently to use the relative number densities (ionization fractions): xm = 
nm/nu - relative abundance of neutral hydrogen, xmi = nuu/nn - the same for ionized hydrogen, a;Hoi = nuci/nuc, xucii = 
nHcii/wHc, a;Hcin = n-Hoiii/nHc - relative abundances of neutral helium, singly and double ionized helium atoms, Xc = ne/nn 

- relative number density of electrons. The ratio of total number densities of helium and hydrogen nuclei we define as /ho = 
WHc/nH, which can be expressed via mass f raction of primordial helium Yp, so that /hb = Yp/4:(1 — Yp) (further we assume 
Yp = 0.24 from ISchramm fc Turner! ilflflSf l'l. This quantities obey evident relationships: Xc = xuu + fneXneii + 2/Hoa;Hciii, 
xm + a^Hii = 1, XHei + XHeii + iHeiii = 1. AH foUowing formulae and relations are presented for these relative number densities 
of atoms, ions and electrons. 

A s it follows from above cited papers and, in particular, from refined numerical calculations by ISeager et al.l il999l 
l20fl(ll at early stages of universe evolution (z > IC*) all hydrogen and helium atoms were ionized completely by thermal 
photons, so, then xmi = 1, xm = 0, ajHoiii = 1, 2;hcI = a^Hcii = and x^ = 1 + 2/ho. This is a consequence of high number 
density of thermal high-energy photons capable to rend electrons from all atoms and ions. In the expanding Universe the 
energy of each photon and temperature of radiation decrease oc a~^, radiation energy density oc a"*, the number and mass 
density of baryons and dark matter oc a~^, where a is the scale factor which is associated with redshift z by simple relation 
a — [z + 1)~^. Already at 2 ~ 8000 thermal photons with energies higher than potentials of ionization of Hell from ground 
and second levels reside in the short-wave tail of Planck function and their number density becomes too low to keep all helium 
in the ionization stage of Helll. It begins to recombine and at 2: ~ 7000 the Hell ions appear. At this moment the time-scale 
of Thomson scattering (tx ^ 3mec(l -I- a;o -I- /Hc)/(8crTaRri|a;o), the hydrogen recombination time-scale [tm ^ l/ricam) and 
the helium one (tnoii — 1/ncaHoii) appear to be essentially lower in comparison with the time-scale of Universe expansion 
(^Hubble — 2/3-ffo(l + z)^'''^)- Thus, matter temperature (electronic and ionic one) Tm equals to CMB temperature, Tr. 
Recombination of HeH occurs in the conditions of local thermodynamic equilibrium (LTE). In the expressions for time-scales 
rrie denotes the mass of electron, c is the light speed, (tt is the effective cross-section of Thomson scattering, br is the radiation 
constant, ai is the effective coefficients of recombination to the ground states of hydrogen atoms HI and singly ionized helium 
atoms, Hell. An ionization fraction of helium, 2;hcIii, is described by Saha equation: 

a;Heii h^riH 

Whereas at this epoch both hydrogen and helium atoms are completely ionized (a;Hi = 0, xmi = 1, a;Hci = 0), a;Hcii = 1 — 2;Heiii 
that means Xc = 1 + ,fuc{l + a;Hein) and above equation can be easily solved for Xc. It proves that already at 2 ~ 5000 all 
helium atoms become singly ionized. This holds up to 2 ~ 3500 when Hel begins to recombine. At this stage fnubbic : : 
tm ■■ tHci ^ 1 : 0.0000003 : 0.0003 : 0.001 and conditions are close to LTE. Until the part of Hel constitutes less than 1% of 
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total helium content the metastable 2s level plays insignificant role in deviation of radiative recombination rate of Hel from 
LTE one and ionized fraction xmi is described yet enough accurately by Saha equation 

Now XHeiii = and a;Hoi = 1 — a;Hoii. To run an accurate calculation of a;Hoii we have to obtain the exact value of = a^nii + 
/He^^Heii. Approximation xnu = 1 is already too rough since 0.1%-decreasing of nmi due to hydrogen recombination results in 
comparable change of rio from Hel recombination because of the prevalent content of hydrogen {fn = nn/iriH + whc) = 0.921). 
So, hydrogen recombination must be taken into account too. At this stage it is described by Saha equation: 

The system of these two equations can be reduced to cubic algebraic equation for Xc which has one real root: 



2^-yl/3cos {a/3)~B/3, (4) 



where B = Rm + Rud, Rud and Rm are right-hand parts of equations lEJ and cos a = 

D = 7?HiiiHoi--RHi-/Hci?Hci, C = 2B^/27-BD/3-E, E = --Rhi-Rhoi(1 - /ho). I have complemented the code RECFAST 
by this solution to achieve correct solution of task formulated above on number density perturbations of ions. However, it 
does not affect the results of calculations for x's noticeably. 

However, such simple description of joint hydrogen- helium recombination loses accuracy very soon, at z ~ 2800. Primarily 
the recombination becomes unequilibrium for Hel because particularity of radiative atomic processes in the expanding cooling 
Universe. Both recombination to the ground state and photoionization from it can be omitted because any recombination 
directly to the ground state will emit a photon with energy greater than potential of ionization and it will immediately ionize 
a neighbouring atom. So, the case B recombination takes place for HeH as well as for HH. Excited atoms of n ^ 2 states are 
ionized by photons of lower energy, belonging to continua of the second and following series, where the number density of 
photons is larger than one of the basic series. Therefore the entire recomb ination proce s s of H el slows down. The cool radiation 
field which is very strong causes the "bottleneck" effect (for details see ISeager et al.| i2nf)ft ll creating the overpopulation of 
excited states relatively to the Boltzmann distribution. So, the Saha equation does not describe adequately the recombination 
more and equations of detailed balance must be utilize d. In the effective th ree-level atom approximation they lead to a single 
differential equation for the ionization fraction of HeH l|Seager et alJll99fll) : 



dxneii _ a^HoiiXennQHcii - /3hci(1 - a:Hoii)e '"^hcI21s/''^"^ 1 + A'HciAHenH(l - 3:Heii)e ^'^v^/^'^'^ 

^ H{z)il + z) 1-f i^Hci(AHc + /3Hci)nH(l-:EHcii)e-'"'^»''"^'" ' 

where 




(6) 



is effective recombination coefficient (m^s~^) of Hel ijHummer fc Storevl llQQj^. /3hcI - photoionization coefficient, KhcI = 
Aggj2ip/[87riii'(2)] - factor which takes into account the cosmological redshifting of Hel 2^p — l^s photons. The values for 
the rest of parameters included into the right-hand part of equation 10 are listed in Table lAll Like in the previous case 
the rate equation for Hel recombination ought to be integrated jointly with equation for HI recombination. Until nni is still 
lesser than 1% of riH the hydrogen ionization fraction, xhii, it can be calculated enough accurate using the Saha equation 
(1600 ^ 2 ^ 2800). Optical depth for Lya emission increases with growth of number of hydrogen atoms in the ground 
state. Diffuse Lya photons, two-photon absorption and collisions as well as cascade recombination from upper levels result in 
overpopulation of the first excited level. Instantaneous spontaneous transition 2p— Is originates Lya photon which is reabsorbed 
immediately by neighbouring HI atoms in the ground state and distribution of levels populations remains unchanged (case 
B recombination). The metastable state 2s of HI is very important for the recombination kinetics because at these redshifts 
the probability of two-photon 2s — Is transition is much smaller than one of photoionization. Therefore, hydrogen atoms 
can be ionized from the first excited state by photons of Balmer continuum, the number density in which exceeds the one 
of Lyman continuum. Jointly with th e overpopulat i on of upper levels relative to a Boltzmann distribution in the strong cool 
radiation field (the "bottleneck" effect ISeager et all i2nf)(l l') it leads to non-equilibrium kinetics of recombinatio n. In this case 
the equation of detailed balance must be used to find hydrogen ionization fraction l Peeble&.1968. : .Seager et allllflflflfl : 

rfa^Hii _ Zoa:HiinHQH - /?h(1 - a:Hn)e~'"^"^°'^'°^"' 1 + A:HAHri-H(l - a^Hii) , 
~dr^ H(z){l + z) i+Kii(AH + l3H)nK{l-xmi)' ^' 

where 

an = F • 10-"atV(l + ct'') m^s"' (8) 
hydrogen recombination coefficient l|Peauignot et al .Il99lh . t = Tm/W^, Kn = X^2p/ [8-^ H{z)] - factor which take into account 
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Figure 1. Hydrogen and helium recombination in ACDM model {Qt = 0.05, ficDM = 0.3, Qa = 0.65, h = 0.65) (top panel). The 
dependences of radiation temperature Tr (K), matter temperature Tm (K), total hydrogen number density nu (m"'^), number density 
of electrons rio (m~^) and Thomson scattering optical depth ''" = /„ caTricH^^dlna on scale factor a = (z + 1)~^ (bottom panel). 



cosmological redshifting of Lya photons. Values for rest of parameters are presented in Table IXTl Photoionization coefficients 
in (O and (O are calculated via the case B recombination coefficients in the following way: 

P = a(27rmefcrn,//i^)=*/'e-''''2--i=/'=^». (9) 

The temperature of electrons and ions, Tm, practically coincides with radiation temperature, Tr, in the range before 
z ~ 800, since until this moment the time-scale of Thomson scattering remains essentially lower than the time-scale of 
Universe expansion, tT/tnubbic < 10~^. Therefore, until this moment the rate of temperature decreasing is governed by 
adiabatic cooling of radiation (7 — 4/3) caused by Universe expansion: 

— = — ^ (10) 

dz (1 + 2) ^ ' 

After recombination, at z < 800, adiabatic cooling of ideal gas (7 = 5/3) begins to dominate over the heating caused 
by Compton effect which is a main process of energy transfer between electrons and photons. Cooling of plasma via free- 
free, free-bound and bound-bound transitions and collisional ionization as well as heating via photoionization and collisional 
recombination gives insignificant contribution into the rate of temperature change, it does not exceed the 0.01% of main 
processes - adiabatic cooling and heating by Compton effect (Seaeer et al. 2000). So, at this epoch the following equation for 
the rate of temperature change proves to be enough accurate l|Wevmann..l96fi,: .Peebles 1968: Seager et al. 1999/1 : 

dz 3H{z)(l + z)m,C 1 + /ho + (1 + 2) ' ^ ' 

The Table IXTl lists the values of all atomic constants and coefficients of approximation formulae seen in equations lfT|l- Hll|l . 
The results of calculations of ionization history in the ACDM model of the Universe, performed on the base of equations 
E^nd RECFAST code complemented by solution JIJ, are presented in Fig[2 There are also shown the dependences 
of radiation temperature Tr, matter temperature Tm, number density of hydrogen nuclei nn, number density of electrons ric 
and optical depth r due to Thomson scattering by electrons t{z) = caTne{z)H~^{z){z + l)~^dz on redshift z (top abscissa 
axis in figure). In consequence of rapid expansion of the Universe and non-equilibrium kinetics of recombination- ionisation 
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processes the hydrogen and helium recombinations do not finish with absolutely neutral medium and some ionized fractions 
persist at low redshifts that is called the residual ionization. The calculations show that at z = 200 ~ xhu ~ 6.7 ■ 10~* 
(FigQ, a;Hoii ~ 9.2 • 10"^" and at z — Xc ~ xmi = 4.1 • 10"'', a;Hon ~ 8.1 ■ 10~^° (an emergence of hydrogen molecules H2, 
and negative hydrogen ion s, H~, does not cha nge these values essentially since their number densities are of several orders 
lower than free electrons ones J Sea.er et aL l Bi . Residual values of ions fractions decrease with increasing of total baryon 
content. In fact, if fib = 0.06 then at 2 = 200 Xe ~ xmi = 6.3- 10~*, xueii = 2.0- 10"^". Stronger dependence of residual values 
of helium ionized fraction (Hell) on baryon content when compared to the hydrogen one (HII) is explained by different atomic 
structures of HI and Hel atoms and conditions under which these recombinations occur: helium due to its higher potential of 
ionization starts to recombine earlier when number densities of plasma particles are higher and its recombination proceed in 
the conditions of high number density of free electrons caused by the complete hydrogen ionization. 



2 PERTURBATIONS OF NUMBER DENSITIES OF IONS AND ELECTRONS 
2.1 Definitions 

Let the values of ionization fractions of HI, HII, Hel, Hell, Helll and electrons averaged over the whole space at fixed 
cosmological time be Xi, where "i" marks each component among them. Let us denote as Xi the local value of relative number 
density of each component in the range of cosmological density perturbation of baryonic matter 5b = Spb/pt <S 1, were pb is its 
matter density. Its deviation from mean value we mark as 5xi, so that Xi = Xi + 5xi and 5xi is called the perturbation of relative 
number density of i-th component. Relative perturbations of relative number densities of ions and free electrons we define 
as Ai = 5xi/xi. It is obvious that Ac — Sne/rie — Snn/nu, Ahii = Snmi/nuu — Snn/riH, AhcII = SnHcu/nucU — Snuc/riHc, 
AhcIii = Snncui/nncui — Snuc/nuc- We suppose that primordial chemical composition of baryon matter is uniform (/hc is 
constant), so Snu/nii = SnHe/nne = 5b and 

A^^5^~5b, (12) 

where Si = Sni/rii = A,; + Sb is relative number density perturbation of i-th component. It must be noted, that in expanding 
universe the recombination does not end with the completely neutral hydrogen or helium, but with residual ionization (see 
last paragraph of previous section). Therefore, none of values rii does reach zero and ambiguity of "0/0" -type in Si does not 
appear. So, A^'s do not take diverging values ever as it seen from 11211 . Numerical results presented below prove that. 

Therefore, Ai is difference of two relative perturbations: of number density of "«"-th component and density of all baryonic 
matter. Since Si and Sb are scalar functions of four coordinates in some gauge, so under the gauge transformations which do 
not change the cosmological background characteristics (mean CMB temperature, isotropic Hubble expansion etc.) each from 
them is t ransfo r med by adding the same expressi o n fro m component of transformation of time coordinate (see for details 
iBardeenl il98(t l: licodama fc Sa,sa,ki| il984 l: IPurreJ i2f)nir i'l. Whereas in 11211 they appear with opposite signs then Ai's keep 
unchanged under such transformations, so they are gauge-invariant quantities. 

If ionization degree does not change with time, that is valid when hydrogen or helium are entirely ionized, then Si — Sb 
and Ai = 0. If the photorecombination and photoionization rates as well as ionization degree of some component change, then 
Si and ^i, can vary with different rates because the altering of Sb is driven by gravitation and stress of baryon-photon plasma, 
while Si is additionally infiuenced by kinetics of ionization-recombination processes. Therefore, Ai is measure of deviation of 
relative number density perturbations of "i"-th component from relative density perturbation of total baryon component as a 
result of changing of recombination and ionization rates within cosmological density perturbation. 

At enough early stage of evolution of the Universe for adiabatic density perturbations of scales larger than the horizon 
Sm, ~ Sb and relative perturbations of radiation energy density 5r = 5eR/eR = 45i,/3. Since ea = uTr, then = STr/Tr = 
l/3Sb. For isothermal perturbations Str = 0. 



2.2 Equations 

Long before and much after recombination the local baryon mass density fiuctuations most probably lead to corresponding 
perturbations of number densities of ions and electrons. Si ^ Sb- But at the recombination epoch because of dependence of 
ionization-recombination process rates on density and temperature of baryon matter the distribution of atoms over ionization 
stages within those perturbations will somewhat depart from background one and Ai / is expected to be true. We study 
the cosmological perturbations of small amplitudes, so the ratio tnubbie : : tei : tuci remains practically the same as 
for background. It means that within cosmological perturbations the same equations ini- Hll|l are applicable and connection 
between the perturbations of ion number density and cosmological perturbations of density and temperature may be obtained 
by variation of those equations. 

Varying the variables wh, T™, a;Hcii, a;HGin and Xg. in the equation JJl we will obtain: 



a;o(l — a;Hoin) 

ZiHcF 



a^o + (1 — a;Hoin)a;Hcin/Ho 



(I 



^ XHcII 



AhcII = — AhcIII, Ac — /hcAhcIII, (13) 



Here it is assumed that xm = a^Hci = at z > 3500 and relative perturbations of rest of components have vanished too. 
One can see that relative perturbations of relative number densities of ions Helll is the linear combination of cosmological 
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perturbations of temperature and mass density of baryon matter. Within an adiabatic perturbation AhcIii has the same 
sign as temperature fluctuation and opposite to mass density one. The values of and iHeiii are calculated from The 
asymptotical behaviour of AhcIii follows from 11811 : at z > 7000 when aiHciii 1 (all helium atoms become double ionized) 
AhcIii — > (feoiii ~ 5b)- On the other hand, a.t z < 5000 when a;Hoin — > AhcIii St„^ and increases with temperature 

decreasing. It is obvious that the second asymptotics is not physically correct. Indeed, such monotonous increasing of AhoIii 
is caused by vanishing of 2;hoIii (see FigQ and does not describe the real number density perturbation of Helll ions. That is 
why in Fig|2| in the range of 5000 < z < 7000 the absolute values of perturbations of number density Sxneiii and Sxuaii are 
presented. 

At 3500 < z < 5000 both hydrogen and helium are entirely ionized: xmi = a^Hcii ~ 1, xm = xuei = a;Heiii = 0. So, at this 
period the amplitudes of all relative perturbations equal to zero. With subsequent decreasing of temperature the Hel atoms 
and afterwards HI ones begin to recombine. The kinetics of their recombination is described by Saha equations Q and 10. 
Variation of these equations gives the following expressions for relative perturbations of relative number density of helium 
AhgII, hydrogen Ahii and free electrons Ac: 



Ahoii = (1 — a;Hcii) 



(1 - a;Hii)a;Hii - fff ) 5t.„ + x, [(§ + fff ) 5t„, - 5b\ 

(1 — 3;Hii)a;Hii + (1 — a;Hcii)a::Hcii,/"Hc — Xa 



(14) 



(1 - xmi)xa [( 



•1 4- XHI 
2 ~^ fcT™ 



)5i 



(1 — a;Hii)a;Hii + (1 — a;Hcii)2:Hcii,/Hc - x^ ' 



(15) 



Ae = 



(1 - xmi)xmi^^^ (1 + (1 - SHciOaJHoii/aJc) 
(1 — a::Hn)a;Hn + (1 — 2;hoIi)2;hcIi/hc — Xc 



+ 



(1 - 2;Hoii)a;HciiffiJ (/hc - (1 - xmi)xmi/xc 
(1 — xmi)xmi + (1 — a;Hcii)a;HGii/Ho — Xc 
(1 — xmi)xmi + (1 ^ a^HciQa^Hcii/Hc 
(1 — xmi)xmi + (1 — a::Hcii)a;Hoii/Ho — Xc 



+ 



(16) 



Their asymptotic behaviour for a^Hcii 1 and xmi — > 1 agrees with our anticipations and appears to be the same as for 
AhcIii case: AhcII, Ahii and Ac — » 0. Other asymptotics for ihcII and 2;hii have not physical sense, since for 
3;hcII ^ 0.99 and xmi ^ 0.99 it is necessary to use the non-equilibrium rate equations and energy balance ll^- llllll . In this 
case the differential equations for relative perturbations Ahii, AhcII and (5t„, are obtained by variation of isll- llllll : 



dAi- 



a::Hcii- 



(l + iCHciAHcnH(l — a;Heii)e 



a;Hoiia;GnHQHcii 



dz 



H{z)(l + z){l + JChci(Aho + /3hci) nH(l — iHciOe 



-hl/ps/feTn, 



Ac + AhcII + Sb + 



CfHcI 



^1 + KuclAHerlHil — 2:Hcii)e 



)/3hoi(1 



a;Hcii)e 



+ 



H{z){l + z){l + A'hoi(Ahc + /3Hci)nH(l - XHcii)e-'^"''=/'=^-) 
dxHcii A'HciAHonH(l — 2;Hcii)e"'"'''='^'°^'" 



5f3B 



a^HcII . 'I'^HoKls r 

AhcII H r7= Ot„ 



dz 



1 + 7('HciAHonH(l — a;Hcii)e 



-hvps/kT^ 



xmi . , hUps J 

AhcII + , ^ OTii, 



1 - a;Hii 



1 — 2;hcii 

da;HGii 



A 



HcII- 



dz 



dxHcii j^Hci(AHc + /3hci)wh(1 - a:Heii)e" 
dz 1 + _K'hci(Ahc + /3Hcl)nH(l - 2;Heii)e 



-hVps/kT^ 



-hl/pa/feTn, 



a^Hii . hVp, 

-AHcII + , m "Tm + 



/3hc 



1 - a;Hii 



Aho + /3hcI /3hcI 



+ 



(17) 



rfAnil _ (1 + -gHlAHWH(l — XHIl)) gHIiacnHaHII [a 4- A 4- A 4- 

dz - H{z)(l + z)(l + KuiiAn + Pm)nHil - xmi)) I ° ' 
(1 + A:HiAHnH(l ~ a;Hii)) Pmjl ~ aHii)e~''""^-/'' '" 
H{z){l + z){l + A'hi(Ah + (3m)nu{l - xmi)) 



5a 



'HI 



QHI 



+ 



'5 13, 



<HI 



5h- 



^ dxHii -K'HiAHrtH(l - a^Hii, 
dz 1 + A:HiAHnH(l - a^Hii 

dxHii -K"hi(Ah + /3Hi)nH(l — xhii) 
dz 1 + A"hi(Ah + /3Hi)nH(l - a;Hii 



xmi 



- xmi 



Sb ~ 



-A 



HII 



-A 



/3hi 1 
dx 



Ahii H — — — di 



- Xmi 



HII- 



'HII 



dz 



a^Hii 



a;Hii 



-Ahii + 



/3hi <?/3h 
Ah + /3hi /3h: 



kTrr. 



+ 



(18) 



dSj 



where 



dz 

27m c- 

— dq 

1 + Z 

dxucu 



3H{z){l + z)meC 1 + /ho + Xc 



4(r,n - Tk)5t,^ + 



dz 



dxi 



~3J^' J f '^" ^^'^ '^'dT '^^'^'^ same as the right-hand side of lO, Q and Ijlljl respectively. The variations of 
recombination coefficients are expressed via temperature perturbations by the equations: 



1 + hi 



1 + /hc + X. 



■{Tm — rR)Ac + Tm^Ti-n — Tr^Tr 



+ 



(19) 



Sam 
am 



d-c-t'' 



Perturbations of ionization fractions 7 




"Hoi 2 \^ 1 + ^Tm/T2 1 + ^Tm/Tl 



) 



Then the variations of photoionization rates can be calculated in the following way: 

5(3i Soi 3 hiyi2s J. 




So, equations Ill7ll - H19|l constitute the system of three ordinary linear differential equations of first order for relative 
perturbations of ions and electrons relative number densities and matter temperature, which can be solved using the publicly 
available code DVERK^. The initial data for them are equilibrium values of relative perturbations of ion relative number 
density at the moment when xmi > 0.99 and a;Hcii > 0.99 calculated by H14II - II15II . 

Let us use the equations lll,'ill -H19 |l to analyse the evolution of relative ion density perturbations and temperature of 
baryonic matter. Sin ce all these equatio ns have solutions in unperturbed problem therefore it seems naturally to supplement 
the code RECFAST iSeager et alJll999fl with block for calculation of perturbations of ionized fractions. The complemented 
code drecfast.f ^ is used further in our analysis of perturbations of ion number density and matter temperature. 



To estimate the magnitude of probable effect it is useful first to consider the stationary adiabatic mass density perturbation 
of baryon matter with some amplitude. Its spatial shape does not matter for our analysis. One may assume for simplicity that 
it is homogeneous in some region of space. 

The results of calculations of ion number density relative perturbations Ill3ll -lll9 | l caused by adiabatic positive matter 
density perturbation (overdensity) with Sb = 10"* and 5t^ = (5i,/3 as well as isothermal {5b — 10"'*, (Jtj, = 0) and "thermal" 
{5b — 0, (5th = ^ • 10~*) ones are presented in Fig[2]for the range of redshifts 200 < z < 10000. One can see that the 
appreciable deviations of relative perturbations of ion number density Ahii, AhoII and Ac from perturbation of total baryon 
mass density ^6 arise in cosmological recombination epoch (800 < z < 1500) as result of amplification of photorecombination 
and photoionization rates. So, the amplitudes of relative perturbations of number densities of protons and electrons are of ~ 5 
times higher than the amplitude of total baryon number density perturbation. For helium ions Hell such ratio is even higher, 
~ 18. However, such large amplitude practically does not contribute to the number density perturbations of free electrons: in 
the Fig[2| Ahii and Ac are superposed in the range of 200 < z < 1500. This is because of too low number density of helium 
ions. Indeed, at z ~ 1200, when xmi — 0.5, iHeii < 0.0003. Also one can see, that relative temperature perturbations are equal 
to the initial ones because of dominant CMB radiation in the energy balance of recombinational plasma. So, the fiuctuations 
in ionization-recombination rates caused by cosmological adiabatic perturbations do not result into additional noticeable local 
fiuctuations of CMB and electron-ion temperatures. 

The figures in the middle and bottom panels illustrate the role of the initial mass density perturbation and initial 
temperature perturbation in arising of perturbations of number densities of ions and electrons. 

For the adiabatic negative (underdensity) perturbations as well as for the isothermal underdensity and cold "thermal" 
ones the pictures are symmetric as follows from equations Ill3ll -lll9 |l . 



3 PERTURBATIONS OF NUMBER DENSITY OF ELECTRONS IN ACDM MODEL 

In previous section we have analysed the perturbations of number density of ions in the framework of "toy" model of stationary 
perturbations. In reality, the amplitudes of adiabatic cosmological perturbations evolve due to gravitational attraction and 
repulsio n by stres s in ba r yon-photon plasma in po te ntial wells created by the d ensity perturbations of dark ma t ter, se e for 
example |Barik£3 fliil): iKodama fc Sasakil flisl: iMa fc RertschinP^eH Tfflfl^: IHu fc Sumvamal ffflfl^: IPurreH EoB) and 
citing therein. When the scale of perturbation becomes substantially smaller than scale of acoustic horizon (Jeans scale), 
then adiabatic perturbations in the baryon-photon plasma start to oscillate like the standing acoustic waves. In consequence 
of recombination the Jeans scale drops and the previously oscillating amplitudes of perturbations in baryon component 
start to monotonously increase mainly under infiuence of gravitational attraction of dark matter density perturbations. The 
amplitudes of perturbations with scale larger than acoustic horizon at recombinatio n epoch increased 5b oc t^^^ in radiation- 
domi nated epoch and 5b o c t^^^ a fter recombination in dust-like Universe. In papers jRardeenl ll98(il: iKodama fc Sa "^ll 9841: 
iMa fc Bertschinger 1995t IHu fc Suerivama 1995; Durrer 20q3) one can find analytical solutions of relevant equations for 
evolution of relative density perturbations in simplified cases of single component mediums as well as the numerical solutions 
for realistic multi-component Universe. 

I shall use here the numerical approach bv lMa fc BertschingeJ 1199^ ) and their package of FORTRAN programs COS- 
MICS'^ in order to calculate the amplitude of baryon density perturbations in multi-component medium for synchronous gauge. 
The evolution of amplitudes of adiabatic density perturbations for each component of Hot plus Cold Dark Matter (HCDM) 



2.3 Results 



^ It is created by T.E. Hull, W.H.Enright, K.R. Jakson in 1976 and is available at site |http:/ /www.cs.toronto.edu/NA/dverk.f.gz| 

^ available at http://astro.franko.lviv.ua/~novos/ 

^ It is available at site http: / / arcturus.mit.edu/cosmics/[ 
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Figure 2. The relative number density perturbations of ions of helium, hydrogen and free electrons caused by adiabatic positive 
(overdensity) matter density initial perturbation (top panel), by isothermal perturbation (middle panel) and by hot "thermal" initial 
fluctuation (bottom panel). 



and ACDM models are shown in the FigEl The scale of these perturbations in Fourier space is fc = O.lMpc"^, that actually 
equals to the horizon distance at the epoch of radiation-matter equality and amounts approximately 1/4 of the horizon size 
at the decoupling epoch. So, starting from radiation-matter equality epoch (a^g ~ 2 ■ lO'^Slm^^) the amplitudes of density 
perturbations in baryon-photon plasma oscillate acoustically till the recombination {arse ^ 10~^). The detailed analysis of 
the ev olution of density perturbations of different scales in the multi-component Universe one can find in lMa fc Bertsching eJ 

In order to analyse the perturbations of number densities of ions and electrons caused by cosmological baryon density 
perturbations the code drecfast.f was complemented by the COSMICS' code linger_syn.f as subroutine, so, that amplitude of 
baryon mass density perturbation 5t and radiation temperature (5th = 5r/4 were precalculated at each step of integration of 
the equations system Ill3ll - lll9|l . Results of joint calculations the ion number density and mass density perturbations of different 
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Figure 3. Evolution of amplitudes of adiabatic density perturbations for different components of HCDM {Cli, = 0.05, f^cDM = 0.75, 
Qi, = 0.2, h = 0.65) and ACDM (Ui, = 0.05, Ocdm = 0.3, S7a = 0.65, h = 0.65) multi-component models in synchronous gauge as 
calculated by COSMICS code (arbitrary normalization). The lines represent the relative density perturbation of baryons Si, (solid line), 
thermal electromagnetic radiation <5a (dash-three dotted line), cold dark matter <5c (dash-dotted line), massless (long-dashed line) and 
massive (short-dashed line) neutrino. The scale of perturbations in Fourier space is fc = O.lMpc"^. 



scales in ACDM model are presented in FigQ] In the left-hand column the evolution of relative mass density perturbations 
5b, 5k and 5c for the baryon, photon and CDM components (adiabatic initial conditions) are shown. Five wave numbers 
are plotted: k — 0.01, 0.025, 0.05, 0.075, 0.1 Mpc~^ (from top to bottom). In the right-hand column the perturbations of 
number density of electrons against baryon mass density ones are shown for the same wave numbers. In each figure the visibility 
function dr /dze~'^ is shown by dotted line in order to mark the position of last scattering surface. The peak of this function 
is found aX z — 1088 and, in fact, denotes moment of decoupling of photons and baryons. The particle horizon at this moment 
equals rjdec — 278Mpc {kdec — 0.023Mpc~^), sound horizon (or Jeans scale) is A^^;, ~ 160Mpc {k^ec — 0.039Mpc~^). Therefore, 
the wave numbers of particle and sound horizons at decoupling epoch fall within the range of plotted k and behaviour of 
baryon mass density and free electron number density perturbations of scales larger, comparable and lower in comparison 
with horizons are revealed. For modes with k < O.OlMpc"^ the relative amplitude, shape and position of peak of electron 
number density perturbation are the same as for k = O.OlMpc"^ mode. For them the peak of 5c is situated at 900 z < 910 
and ratio of amplitudes 5c/5i, ~ 4.5. For modes with k > O.OlMpc"^ the relative amplitude, shape and position of peak of 
electron number density perturbation vary: when a wave number increases starting from 0.01 to k^^c the peak position 5^ 
shifts to the position of visibility function peak with approximately the same ratio of 5c /5b ~ 4.5. For modes with k > k^cc 
the value of 5c radically decreases and makes several oscillation around value of 5b depending on phase of (5r oscillation at 
decoupling epoch. So, ion number density fluctuations at decoupling epoch caused by recombination kinetics in the range 
of adiabatic density perturbations are most prominent for the large-scale perturbations with k ^ k^cc- At z < 800, when 
ionization become residual, 5c becomes lesser than 5b independently of scale of baryon density perturbation. The difference 
acquires the value ~ 20% at z 200. It is caused by increasing of 5b on all scales under gravitational potential of dark matter 
density perturbations, so, the probability of recombination is slightly higher in such regions. It is similar to the residual 
ionization and its inverse dependence on baryon density. 

Finally it should be noted that presented results do not depend essentially on exact values of cosmological parameters of 
non-exotic models and are practically the same for parameters of WMAP concordance model. 




Figure 4. Left-hand column: evolution of relative density perturbations Sf,, i5r, and 5c for the baryon (solid line), photon (dash-three 
dotted line) and CDM (dash-dotted line) components with adiabatic initial conditions and free normalization. Right-hand column: 
electron number density perturbations against baryon mass density ones. The amplitudes of all perturbations are in dimensionless units 
10~®. The visibility function dT/dze~'^ (multiplied by 0.018 for convenience, dotted line) denotes the position of last scattering surface. 
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4 Ae and corrugation of last scattering surface 

The results presented in Fig[2] prove that ion number density fluctuations caused by kinetics of ionization-recombination 
processes in the field of cosmological adiabatic perturbations do not lead to appreciable additional local temperature fiuctu- 
ations. This was expected because of dominating role of the thermal (relic) electromagnetic radiation in the energy balance 
of photon-baryon plasma. However, to calculate of CMB temperature fluctuations and polarization at high accuracy in order 
to achieve better agreement between theory and observations it is necessary to take into account the contributions of such 
electron number density fluctuations into optical depth r caused by Thomson scattering . 

The visibility function dr / dze~^ , where t{z) = f^ caT'nc{z)H^^{z){z + l)^^dz, represents the probability that a photon 
was scattered for the last time within dz of z. So, the main part of CMB photons come to observer from the vicinity of maximum 
of visibility function. The fluctuations of number density of electrons result into faint "corrugation" of last scattering surface of 
thermal relic radiation: the maximum of visibility function will be at somewhat lower redshifts for overdensity perturbations 
and at somewhat higher redshifts for underdensities than for unperturbed region. We can estimate this effect in the following 
way. 

The redshift of peak can be determined from the condition of local extremum giving the equation: 

1 dxs _ 1 dH(z) 2 _ caTUHXeiz) _ 

^~d^~ H{z) dz 7TT^ H(z){z + 1) ^ ' ^ ' 

where ne{z) is unperturbed number density of electrons shown in Fig[2 H{z) = 

Ho ^Q-r{z + 1)''' + Qm{z + 1)^ + Qk{z + 1)"^ + ^ is Hubble constant. Its solution for z gives a position of peak of 

visibility function in unperturbed medium. For ACDM model with the same parameters as for Fig|l| using the numerical 
solution we obtain Zdec — 1088. In the range of perturbation, where fiu — nH(l + 5b) and = Xe{l + Ao), it is expected 
to be Zdec = Zdec + Szdec- The displacement of peak Szdec can be estimated by variation of equation H2()ll and expanding of 
functions Xe{zdec) and nii{zd^c) into Taylor series about Zdec (in linear approximation): 

CaTUHXeiz) dAc 
(Ot + Ac 



Z+1 



H{z)(z + 1)^ " ^' dz 



caTTiuXeiz) ( 1 dH{z) ^ 1 ^^ ^ '^^•^ ( l1^o^_^-|^^^/^ ^ dH[z) 1 dxr, 

z + i) (z + i) I \ {z -\- I) 



H(z)[z + 1) \H{z) dz Xe dz J dz \H{z) dz Xe dz J z + 1 



(21) 



Therefore, such displacement depends on the amplitude of Ac + St = Sc and the gradient of Ac. For arbitrary normalisation 
of amplitude of cosmological perturbations it is conveniently to present the result in units of relative density perturbations 
of baryon-photon plasma components and for large-scale perturbations {k < kd^c) we have: Szdec/ {zdec + 1) « — 0.255r = 
—0.33(56. If we suppose Ao = (ion number density perturbations follow the baryon mass density ones) then Szdecj {zdcc + 
1) — 0.0515r = —0.068(56. So, the amplification of perturbation amplitudes of electron number density by recombination 
processes makes the last scattering surface more "corrugated" in optical depth. Such "corrugation" results into observable 
CMB temperature fluctuations. 

We can make the estimation of this effect. Let us suppose that major part of CMB photons come from thin last scattering 
surface placed at Zdcc from unperturbed medium and idee = Zdec + Szdec from perturbed one. Because temperature of 
the thermal radiation at Zdec and at idee is the same then its observable variation follows from the well-known relation 
To = Tdec/{zdec + 1): {STo/Tq)^^^ — —5zdec/{zdec + 1) « 0.25(5r = 0.33(56 that equals to intrinsic adiabatic CMB temperature 
fluctuation, ((5To/To)^^ — 5r/4 — ^6/3, exactly. It is other way for calculation of adiabatic term of CMB primary anisotropy, 
which, however, requires the accurate precalculation of Ae or 5e different from 5b at last scattering surface. When we suppose 
that &e = 5b {Ae = 0) then (5To/To)^^^ « 0.2 {5To/To)^^ that is incorrect. 

For small-scale perturbations k > kdeo the value of 5zdec/ {zdec + 1) will be smaller and its sign will alternate depending 
on oscillation phases of (5r and 5b at Zdec- The assumption of thin LSS is too rough in this case to make correct estimation. 
The line-of-sight integration must be undertaken here because of fuzziness effect leading to exponential reduction of any 
primary CMB anisotropy. Unfortunately, straightforward substitution of recfast.f subroutine by drecfast.f one in available 
codes (CMBfast, CMBEASY etc.) for calculation of temperature and polarization Ci's does not give correct estimation of 
possible changes to CMB power spectra caused by effect demonstrated here. It looks that collision term connected with 
Thomson scattering in electron density perturbed region and Boltzmann equation for photons must be generalized to take 
this effect into account properly. It will be the matter of the separate paper. 



conclusions 

In the fleld of adiabatic density perturbations the rates of hydrogen and helium ionization-recombination processes slightly 
differ from those in the non-perturbed medium. On the one hand the rate of recombination increases due to somewhat larger 
number density of plasma particles in the overdense region, on the other hand the photoionization is increased too due to 
positive temperature initial perturbation since number density of photons capable to ionize atoms is slightly higher, rising the 
level of ionization. The second effect prevails at the beginning of recombination epoch and competition of effects results into 
the positive additional overdensity of ions and electrons. For initial perturbations with opposite sign (underdense region) all 
results are symmetric. The effect is prominent for large-scale initial perturbations {k ^ kd^c) which never oscillate acoustically. 
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Thus, the maximal amplitude of perturbations of proton and electron relative number density during the recombination epoch 
is by factor of ~ 5 higher than amplitude of Sb, baryon density perturbations. For helium it is ~ 18 times higher, but practically 
does not influence the amplitude of perturbations of electron relative number density because of low level of ionization at 
z < 1200. For initial perturbations of scales inside acoustic horizon at decoupling epoch (fc ^ fc^ec) the deviations of ion 
number density perturbations from baryon total density ones vanish owing to the temperature and density oscillations of 
photon-baryon plasma during recombination. They are small also for isothermal initial perturbations. 

At lower redshifts, when ionization becomes residual and baryon mass density perturbation increases the electron number 
density perturbations 6e becomes lesser than S), independently of scale of baryon density perturbation. The difference acquires 
- 20% at 2 ~ 200. 

Revealed deviations of electron number density perturbations amplitudes from baryon total density ones at cosmological 
recombination epoch do not lead to appreciable additional local temperature fluctuations in matter or thermal radiation. But 
they result in faint optical depth "corrugation" of last scattering surface: Szdec/ {zdec + 1) ~ — 0.33(5;, at scales larger than 
sound horizon. Taking them into account may improve the agreement between theoretical predictions and observable data of 
current WMAP and future PLANCK missions on large-scale CMB anisotropy and polarization. 
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Table Al. Atomic constants and coefficients of approximation formulae 



Constant Value Source Formula 



XHI 2.17871122 ■ 10-i» J Seager et al. (im9) ISJ,CSJ 

XHoi 3.9393393 • 10-1** J Seager et al. fl999l 0, {HJ 

XHoll 8.71869443 ■ 10-i» J Seager et al. 

hun2s 1.63403509 ■ 10-1** J Seager et al. (19991 ffl.lTsl 

hups 3.30301387 ■ 10-1** J Seager et al. (1999) ©,C2l 

f^'^UcU^s 9.64908312 ■ 10-20 J Seager et al. fl999l isl.lTTt 

hu2s-is (HI) 5.4467613 ■ lO-^'' J Peauignot et al. (1991) ^ 

hu2s-.is (Hell) 6.36325429 ■ lO-i^ J m immer fc Stor ey (1994) 
Ah2p 121.567 nm Seager et al. (19991 : Verner fc Ferland (19961 0,CHJ 
AhoI21p 58.4334 nm Seager et al. (19991 : Verner fc Ferland (19961 isl.lTTt 

F 1.14 Seager et al. (19991 

a 4.309 Pennignot et al. (19911 |H) 

b -0.6166 Peauignot et al. (19911 |8j 

c 0.6703 peauignot et al. (iggl) HJ 

d 0.5300 peauignot et al. (19911 |8j 

q 1.80301774 ■ 10-" ^mmer & Storey (19981 IHJ 

p 0.711 jiummer & Storey (19981 |6j 

Ti 1.30016958 ■ 10^ K Rummer & Storey (19981 IHJ 

T2 3K 'Hummer fc Storey (19981 |6j 

Ah 8.22458 c-1 Goldman (19891 0,CHJ 

Aho 51.3 c-1 Drake et al. (19691 isl.lTTt 



APPENDIX A: THE VALUES OF ATOMIC CONSTANTS AND COEFFICIENTS OF 
APPROXIMATION FORMULAE USED IN NUMERICAL CALCULATIONS 



